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Entropy Crisis, Defects and the Role of Competition in Monatomic Glass Formers 
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Abstract 

The demonstration here of an entropy crisis in monatomic glass formers along with previously known results finally establishes 
that the entropy crisis is ubiqutous in all supercooled liquids. We show that interactions that compete with crystal order weakens 
the ability to form glass. We also study the defects in the ideal glass; they are very different from those in the crystal. 



A common feature of almost all materials is that they 
can become glassy, i.e. an amorphous solid, representing 
a metastable state The glass is obtained by cool- 
ing the supercooled liquid (SCL), and need not obey the 
third law of thermodynamics, which is only valid for equi- 
librium states 0] . The true equilibrium state corresponds 
to a crystalline state. The configurational entropy of SCL 
exhibits a rapid drop below the melting temperature Tm 
and eventually vanishes at some temperature T — 
Tk < Tm. It will become negative if extrapolated to lower 
temperatures. Since a negative entropy is unphysical, an 
ideal glass transition must intervene to avoid this entropy 
crisis at Tk. In the limit of zero cooling rate (not acces- 
sible in experiments or simulations, but accessible in a 
theoretical setup) in the metastable region, metastable 
states become stationary, and can be described by " equi- 
librium" statistical mechanics by restricting the allowed 
microstates to be disordered; microstates leading to the 
crystalline state are not allowed. In experiments, the 
extrapolated configurational entropy of many glassy ma- 
terials at T = [3| is found to have a non-zero value 
depending on the rate of cooling. This does not rule out 
the possibility that the entropy of the hypothetical " sta- 
tionary" glass vanishes in the limit of zero cooling rate 
at a non-zero temperature at a positive temperature. 

Demonstrating an entropy crisis for supercooled liq- 
uids in a restricted formalism has been one of the most 
challenging problems in theoretical physics. An entropy 
crisis in long polymers was theoretically demonstrated 
almost fifty years ago by Gibbs and DiMarzio ^ to sup- 
port the entropy crisis as a fundamental principle under- 
lying glass transitions in long polymers. Their work was 
later severely criticized by Gujrati and Goldstein [5] for 
its poor approximation, and doubt was cast on whether 
the entropy crisis in long polymers was genuine. The 
situation changed when in a recent work, Gujrati and 
Corsi [6] established the existence of an entropy crisis in 
long polymers by using a highly reliable approximate ap- 
proach. This was very important as the idea of Gibbs 
and Di Marzio has been pivotal in shaping our thinking 
about the ideal glass transition. Recently, we have also 
succeeded in demonstrating the entropy crisis in a dimer 
model 7] containing anisotropic interactions. However, 
to establish the entropy crisis as a fundamental principle 
underliving the glass transition in all supercooled liquids, 
we need to demonstrate the crisis in simple isotropic flu- 
ids containing monatomic particles. So far, this has not 



been feasible in any theoretical approach and is one of the 
outstanding theoretical physics problems for a complete 
understanding of glass transition. 

Here, we establish an entropy crisis in monatomic 
systems, thereby finally succeeding in establishing the 
entropy crisis as a fundamental mechanism driving 
glass transitions in all supercooled liquids. Study of 
monatomic glass formers will allow us to obtain a better 
understanding of the glassy structure (defects therein) 
whose accurate representation remains still challenging. 
The slow relaxation jlj in SCLs is similar to that observed 
in ordinary spin glasses Q. Therefore, it is not surprising 
that ordinary monatomic glasses have been traditionally 
modeled as spin glasses, whose important features are 
their geometrical frustration and competition. It is com- 
monly believed that the competition and frustration play 
an important role in promoting the glassy behavior. We 
will also check this hypothesis in this work. 

It should be noted that frustrated antiferromagnets 
(AF) and spin glasses do not usually possess long range 
order at low temperatures because of a highly degener- 
ate ground state |3] and their glassy behavior is brought 
about by the presence of frustration or quenched impuri- 
ties and is somewhat well understood. In contrast, super- 
cooled liquids require a unique ground state, the crystal. 
This distinguishes the glassy behavior in supercooled liq- 
uids and requires considering an unfrustrated AF model 
as a paradigm of simple fluids or alloys. We consider a 
pure (no frustration or quenched impurities) AF Ising 
model, which possesses a unique ordered state so that 
supercooling can occur. This then results in a glassy 
state. We are not aware of any simple model calculation 
to date to justify glassy states in a pure AF model. We 
also flnd that the competition considered in this work in- 
hibits instead of promoting the glass transition, which is 
a surprising result. 

Model. We introduce the following AF Ising model in 
zero magnetic held on a square or a cubic lattice (lattice 
spacing a) with the interaction energy 
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The flrst sum is over nearest-neighbor spin pairs and the 
second over neighboring spin triplets, which we take to 
be within a square for simplicity. We take = ±1 to de- 
note A, and B particles for an alloy, or the particle and 
void for a fluid. For | J'| < 2 J, we have an AF ordering 
at low temperatures with a sublattice structure: spins of 
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a given orientation are found preferentially on one of tt 
two sublattices. Antiferromagnetically ordered square 
(AFS) with spins alternating, and ferromagnetically o: 
dered squares (FS) with spins the same are the only tw 
square conformations that contribute to the first ten 
in ([T]). We may identify the AF ordered structure as 
crystal [7|. For a fluid, this model represents a stron 
repulsion at a lattice spacing a, and attraction at lattic 
spacings \/2a between particles. For | J'| > 2 J, the A 
ordering is destroyed at low temperatures; S is the sarr 
everywhere. We set J=l to set the temperature sea' 
and only consider | J'| < 2J. It is easy to see that th 
free energy depends on | J'| , not on its sign. In partici 
lar, the ground state energy per spin of the AF ordere 
state is £'0 = — 2J, regardless of J'. In the following, w 
will measure the energy and the free energy with respee 
to the ground state to give the excitation energies. I 
this case, both will vanish at T = 0. The non-zero vak 
of I J' I creates a preference for the product SS'S" in 
square to be of a fixed sign, which then competes wit 
the formation of the crystal in which this product ca 
be of either sign. A positive (negative) J' provides 
preference for 5" = — 1 (+1), so J' can be used to als 
control the abundance of one of the spin states. 

The entropy >S'(r) of the model cannot be negative 
if the state has to occur in Nature or simulations; in- 
deed, neither can ever show any entropy crisis. If the 
metastable state entropy S(T) = at a positive temper- 
ature Tk, then its extension will experience an entropy 
crisis and must be replaced by an ideal glass below Tk, 
the ideal glass transition temperature [lO|. The model 
cannot be solved exactly except in one dimension. It 
is usually studied in the mean-field approximation com- 
monly known as the Bragg- Williams approximation 
adapted for an AF case. However, the approximations 
is known to be very crude. Indeed, Netz and Berker 
[12.] have shown that one of the shortcomings of the ap- 
proximation is that it abandons the hard-spin condition 
— 1. This condition is easily incorporated in exact cal- 
culations on recursive lattices [Tsj and it was discovered 
that such calculations are more reliable than the conven- 
tional mean-field approximations. Therefore, we adopt 
the recursive lattice approach here. 

We consider a homogeneous Husimi cactus in which q 
squares meet at a site. The model is solved exactly on the 
cactus as described in [l^. The cactus can be thought 
as an approximation of a square lattice for q — 2 or a 
cubic lattice for g = 3, so that the exact Husimi cactus 
solution can be thought of as their approximate solution. 

Solution. We follow [l^ and solve the model recur- 
sively. We label sites on the cactus by an index m, which 
increases sequentially outwards from m = at the origin. 
We introduce partial PF's Zm{^) and Zm{i), depending 
on the states of the spin at the m-th cactus level. It rep- 
resents the contribution of the part of the cactus above 
that level to the total PF. We then introduce the ratio 

= Z„(t)/[Z„(T) + (2) 
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FIG. 1: The free energy and entropy for the two FP solu- 
tions. The model shows an entropy crisis and an ideal glass 
transition at Tk. The thin curves represent unphysical states 
(negative entropy) and are replaced by the ideal glass state. 



which satisfies the recursion relation (RR) 

f{x„i+l,Xm+2,v) + f{ym+l,ym+2A/v)' 

where f{x,x',v) = x'^^x'^ /uv'^ + 2x'^x'^y'^v + x'^^y'^v + 
ux'^y'^'^' + 2x^y^y'^ + y'^^y'^ /v with r = q — 1 and where 
u = e^^ , V = e^^'^' ,y = 1 — x,y' = 1 — x' . 

There are two kinds of fix-point (FP) solutions of the 
RR ([3]) that describe the bulk behavior [3]. In the 1- 
cycle solution, the FP solution becomes independent of 
the level m, and is represented by x* . For J' = 0, x* 
is given by x* = 1/2, as can be checked explicitly. For 
J' 0, a;* 7^ 1/2 and has to be obtained numerically. 
This solution exists at all temperatures T > 0; thus, there 
is no spinodal of this solution. This solution describes the 
disordered phase. At T — > 00, all spins are uncorrelated 
so that the density per site 4>-ps = 1/8, and 0afs = 1/16, 
and the entropy is 5* = In 2. As T is reduced, (/)fs de- 
creases, while 0AFS increases. The other FP solution of 
interest is a 2-cycle solution associated with the AF state 
containg AFSs [l^. It alternates between two values x'^, 
and X2 which occur at successive levels. This kind of FP 
solution has also been observed in other systems such as 
semi- fiexible polymers [13] , dimers [T5l| , and stars and 
dendrimers [l6|, and has been thoroughly investigated. 
At T = 0, the 2-cycle solution is given by x\,X2 — 1,0 or 
0, 1 describing the perfect crystal (^fs = O,0afs = 0.5). 
This solution then evolves with T due to excitations and 
describes the crystal at low temperatures. The free en- 
ergyis calculated by the general method due to Gujrati 
[1, [13I [T6| . Whichever solution has the lower free en- 
ergy represents the equilibrium state. The solution with 
the higher free energy, then, represents the metastable 
state, which can only be observed in Nature if its en- 
tropy remains non-negative. The temperature where the 
two solutions have the same free energy is the transition 
temperature, which we denote by Tm. 

Results. The results for g = 2 are presented in Figs. 
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FIG. 2: S-E-T relationship for the two FP solutions. The 
excitations in the two near T = are very different. The ex- 
citations in the 1-cycle state near Tk are strongly interacting 
as opposed to those near T = 0. 



1-3. The free energy Fi and entropy 5"! associated with 
the 1-cycle FP solution are shown by the continuous and 
the long dash curves in Fig. 1. The free energy F2 and 
entropy 5*2 associate with the 2-cycle FP solution are 
shown by the dotted and the dash-dot curves. The en- 
ergy E{T) as a function of T and the entropy S{E) as a 
function of E are shown in Fig. 2. We have set J' = 0.01 
for figures 1, and 2. As said earlier, F and E repre- 
sent the contributions of excitations with respect to the 
ground state energy Eq = —2 J, so that they vanish at 
T = 0, as is clearly seen in Figs. 1, and 2. The transition 
temperature is found to be Tm = 2.7706. We see from 
Fig. 1 that Fi crosses zero and becomes positive below 
T = Tcq ~ 2.200 but again becomes zero (not shown 
here, but we have checked it) as T — > 0. Thus, Fi pos- 
sesses a maximum at an intermediate temperature (see 
point K in Fig. 1) at T = Tk 1.1316, so that the en- 
tropy 5*1 vanishes there. Below Tk, the continuation of 
Fi and 5*1, shown by their thin portions in Figs. 1 and 2, 
continue to satisfy the stability condition (non-negative 
specific heat). Despite this, they cannot represent any 
physical states in the system due to negative entropy and 
have to be discarded as unphysical. Below Tk, we must 
extend the metastable state (described by Fi and 5*1 be- 
tween Tk and Tm) by a glassy phase of a constant free 
energy F = Fq, sec the short dash horizontal line in 
Fig. 1, and S ~ Sq ~ 0- The 1-cycle energy at K is 
EiK = Fq ~ 0.301. The entropy S2 is never negative, 
and the 2-cycle FP solution represents the equilibrium 
crystal below Tm. 

Competition with Crystal Ordering. The behav- 
ior as a function of \J'\ of the transition temperature 
Tm (empty circles), the ideal glass transition tempera- 
ture Tk (filled circles) and their ratio Tm/Tk (triangles) 
are shown in Fig. 3. As said above, | J'| competes with 



the crystal ordering and reduces Tm- One can use the 
inverse ratio Tk/Tm as a measure of the relative ease 
of glass formation: larger this value, easier it is to ob- 
tain the ideal glass as Tk is not too deep relative to Tm. 
What we observe is that Tm/Tk increases with | J'| , with 
Tk approaching zero faster than Tm, so that the ratio 
Tm/Tk continues to increase with | J'| . This implies that 
it becomes harder to obtain the ideal glass as Tk be- 
comes relatively farther away from Tm as \J'\ increases. 
The competition provided by \J'\ weakens not only crys- 
tal ordering but also "weakens" forming the ideal glass. 
Consequently, competition does not enhance the ability 
to undergo ideal glass transition, an interesting result 
which is being explored further [l3| to see if other com- 
petitions behave similarly. 

Analysis of Defects. To understand the difference 
between the disordered liquid and the crystal defects, we 
turn to Fig. 2 and observe that near T = 0, the exci- 
tation energies of both FP solutions are very different, 
even though Ei(0) = EM) — ^- Detailed analysis will 
be pesented elsewhere [l3|. The excitations (defects) in 
the crystal are known to be due to point-like excitations 
caused by the reversal of a single spin which changes the 
free energy by ~ 8J (coordination number 4) with re- 
spect to the ground state; here we assume J' to be small. 
Therefore, this excitation causes the leading term in the 
free energy F2 to be 1 /u^ [l^ and can be treated as non- 
interacting as long as they are small in number. What 
kinds of excitations are deducible from the form of Ti 
near T = 0? To answer this, we consider the simple case 
of J' = 0, which also gives rise to Ti over the entire tem- 
perature range T > 0. Here, we can carry out an analyti- 
cal investigation since x* = 1/2. As we are only interested 
in the excitation energy, we will overlook the unphysical 
nature of the entropy near T = to study the excitation 
spectrum. The excitation energy [=4(</)fs — </'afs)+2] 
and free energy are given by (lu = 3 -f u/2 + l/2u) 

Ei{T) = -{u-l/u)/w + 2,Fi{T) = -(T/2)lnw-t-2, 

where 0fs = 1/4mw, and 0afs = u/Aw. The entropy is 
calculated using S'i(T) = {Ei-Fi)/T. At T = 0, Ei{0) = 

implies that all squares are. Near T = 0, we find that 
£:i(T) ~ 12/it, Ti(T) ~ T(iln2 - 3/u) and S'i(T) ~ 
-i ln2 -f 3/u -\- UP/u. {Si at T = in Fig. 1 is almost 
— i ln2.) What we discover is that the excitations due to 

1 /w-term near T = are not the uncorrelated single spin 
reversal in the background of a perfect crystal. Rather, 
they represent correlated excitations in the form of FS 
in the background of AFSs by turning an AFS into a 
FS. Each FS excitation requires an energy 4J per site. 
This is in accordance with the general discussion above. 
Indeed, the excitation spectrum is given by the expansion 
of Ei{T) in powers of 1/m [l^. These excitations also 
explain why the thin portion of Ei (T) rises more rapidly 
than E2{T) in Fig. 2 so that Ek — 0.3 is appreciably 
higher than the E2{T) ~ 0.1 at Tk due to a lower (^afs — 
0.426 and higher 0fs — 3.58 x lO""*^ in the metastable 
state. In contrast, 0afs — 0.498, and (f)Fs — 7.991 x 
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FIG. 3: The effect of \J'\ , which creates competition with 
the crystal ordering, on Tm, and Tk and their ratio. The 
weakening of crystal ordering also "weakens" the ideal glass 
formation. 



10^^ in the crystal. This is consistent with most defects 
in the crystal being point defects. As there is no non- 
analyticty at Tk in SCL, the excitation spectrum remains 
continuous above Tk where physical states occur. The 
density of sites (j)unc not covered by AFSs and FSs is 
0u„c ~ 0.148 in SCL and (ji^nc ^ 3.684 x 10"^ in the 
crystal. These sites are probably uncorrelated in SCL, 
but this needs to be carefully checked. Nevertheless, the 
above analysis leads us to conclude that the excitations 
at or above Tk in SCL are very different from the point 
defects of the crystal and are the ones that get frozen in 
the glass that is formed at Tk- For T < Tk, we have an 
ideal glass, shown by the horizontal short dash curve in 
Fig. 2, of constant energy iJx and zero entropy. 

The excitation spectra of both solutions over the physi- 
cal range are completely described by their respective en- 
tropies Si{E) (the thick part) and S2{E) shown in Fig. 2. 
For E < i?K, the excitations in the glass cannot change 



since they are frozen at constant energy i?K, but con- 
tinue to change in the crystal. In experiments, the ideal 
glass will never be observed due to time-limitations and 
one would obtain a non-stationary state whose entropy 
^1 {E) must satisfy S\(E) < Si {E) according to the law 
of increase of entropy . In this case, the non-stationary 
glass will have some excitations at low temperatures. 

For the AF case that we consider here, the 1-cycle so- 
lution is found to exist at all temperatures and describes 
the disordered liquid above and its metastable contin- 
uation below the transition temperature. There is no 
singularity in this fix-point solution at the transition. In 
contrast, for the ferromagnetic case (J < 0), the 2-cycle 
FP solution is never stable, and the 1-cycle solution has 
a singularity at the ferromagnetic transition and its en- 
tropy never becomes negative. 

We have already shown elsewhere p^] that the glassy 
state for dimers contains a higher density of voids than 
the corresponding crystal at the same temperature. This 
is also true in the current model. The voids distribute 
themselves in the lattice at equilibrium, and the corre- 
sponding 1-cycle solution gives rise to an excitation spec- 
trum so that Si{E) vanishes at Tk. If it happens that the 
system is quenched, then all we can say is that the corre- 
sponding spectrum Si{E) of the quenched system must 
satisfy the standard condition Si{E) < Si{E). Despite 
this, it is possible that the entropy of the quenched sys- 
tem does not vanish at a positive temperature. There is 
no contradiction. 

To summarize, our model calculation demonstrates 
that monatomic systems also give rise to an ideal glass, 
thereby making the entropy crisis ubiquitous. The glass 
contains correlated defects that are very different from 
those in the crystal. Moreover, the competition does not 
enhance the ability to form a glass. We have only inves- 
tigated a classical model and it would be interesting to 
see if quantum calculations support this picture. 
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